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Изучается возникновение, развитие и последствия акустической неустойчивости в камере ракетного двигателя с твердотоплив- 
ным зарядом. Рассматриваются относительно низкочастотные колебания с периодом, намного превышающим период собствен- 
ных колебаний камеры. Но частота изменения термодинамических параметров находится в широких пределах и может быть 
сравнима с собственной частотой зоны горения. Неустойчивость может приводить к автоколебательному горению, или к хаоти- 
ческому режиму, или к погасанию горения. 


Введение 

Еще в конце 30-х гг. прошлого века разработчи- 
ки системы залпового огня «Катюша» столкнулись 
со странным явлением. Как оказалось, ряд ориги- 
нальных технических решений, который позволил 
бы достичь больших скоростей реактивных снаря- 
дов, не осуществим из-за особенностей горения то- 
плива. Оно, превосходно сгорая на открытом воз- 


духе, при определенных условиях «не желало» го- 
реть в камере, или же горело настолько нестабиль- 
но, что совершенно приводило в негодность двига- 
тель. Как оказалось, наличие окислителя и горюче- 
го в составе топлива не гарантирует равномерное 
протекание химической реакции, необходимо еще, 
чтобы внутренние закономерности механизма го- 
рения оптимально сочетались с конструкционны- 
ми особенностями самого двигателя ракеты. Даль- 
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нейшие теоретические и экспериментальные ис- 
следования, проведенные в основном в СССР и 
США [1-9], позволили понять основные причины 
нестабильного горения в ракетном двигателе твер- 
дого топлива (РДТТ) и сформулировать простые 
критерии для избегания проявления неустойчиво- 
сти. Но многогранность и сложность явления горе- 
ния показывает, что эти критерии имеют весьма уз- 
кую границу применимости [1]. Особенно это на- 
чинает ярко проявляться по мере продвижения в 
высокочастотную область, при горении веществ со 
сложной кинетикой разложения и ее изменчивости 
при различных температурах и давлениях [2]. 

Поэтому дальнейшие усилия должны быть на- 
правлены на раскрытие детального механизма го- 
рения, установление его связи с конструкционны- 
ми особенностями РДТТ при возникновении неу- 
стойчивости. Настоящая работа является продол- 
жением комплекса теоретических исследований 
[3-5], направленных на выявление такой модели 
нестационарного горения, которая содержала бы 
минимум параметров и описывала максимум на- 
блюдаемых эффектов. 

Математическая модель 

нестационарного горения в РДТТ 

При возникновении неустойчивого горения в 
РДТТ, как правило, звуковые колебания давления 
происходят на основной моде. При этом не обяза- 
тельно, чтобы частота звуковой волны совпадала 
или была близка к собственной частоте горящего 
топлива: при нарушении баланса массы продуктов 
горения в камере сгорания колебания могут проис- 
ходить и на частоте, близкой к обратной величине 
характерного времени истечения газа из сопла. В 
таком случае длина волны звука много больше дли- 
ны камеры, все термодинамические и гидродина- 
мические параметры газа можно считать не завися- 
щими от пространственных переменных. Это, ко- 
нечно, не относится к процессам горения, которые 
протекают на масштабах порядка нескольких де- 
сятков мкм. Здесь рассматриваются такие измене- 
ния физических величин, которые характеризуют 
двигатель в целом, и их с удовлетворительной точ- 
ностью можно считать зависящими только от вре- 
мени. Это существенно упрощает теоретический 
анализ. 

Пусть 5 Р - площадь поверхности горения; и - 
скорость горения; і— время. Обычно (наибольшее) 
характерное время протекания нестационарных 
процессов таково, что поверхность горения за этот 
промежуток времени меняется слабо. Поэтому 
можно считать Д»сошІ. Скорость изменения мас- 
сы т=рѴ с газа с плотностью р в камере объема Ѵ с 
определяется разностью поступающего за счет го- 
рения и покидающего камеру через сопло количе- 
ства вещества [1]: 

< ^7 = Рс и3 р- А сР р т’ (!) 

ш 


где р с - плотность топлива; Д, - коэффициент 
истечения; Е т - минимальное (критическое) сече- 
ние сопла. Плотность газа можно определять по 
уравнению состояния идеального газа р=р/І{,Т р . 
Здесь р - давление; К, - газовая постоянная; Т р - 
температура продуктов сгорания. В широких пре- 
делах изменения плотности и давления газа незна- 
чительными колебаниями температуры продуктов 
сгорания можно пренебречь (изотермическое при- 
ближение) [1]. Тогда ур. (1) можно записать как 
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Начальным условием для этого уравнения служит 
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Волну горения твердого ракетного топлива, 
движущуюся в отрицательном направлении коор- 
динаты х, представим следующей моделью [4, 5]: 
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Граничные условия: 
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х -> +оо : ^р- = 0, = 0. 

Ах Ах 

Здесь х з - поверхность разложения топлива; с„ 
к — теплоемкость и коэффициент теплопроводно- 
сти материала топлива с температурой Т с , 
Т с (х=х)= 7]; Г 0 - начальная температура топлива; У, 
7) - массовая концентрация (доля) и коэффициент 
диффузии реагирующего вещества; к 0 - предэкспо- 
ненциальный множитель в законе Аррениуса; Е„ 
Е - эффективные энергии активации химической 
реакции в твердой и газовой фазах; Т- температура 
газа; К — универсальная газовая постоянная; Я - ко- 
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эффициент теплопроводности газа; 0 - суммарный 
тепловой эффект химической реакции в газе; у- по- 
казатель адиабаты; Ь - тепловой эффект разложе- 
ния топлива на газообразные компоненты; ѵ 0 =соп$1. 


В качестве начальных условий к ур. (3) должны 
быть взяты их стационарные решения. Дальней- 
шие упрощения для решения ур. (2), (3) сводятся к 
переходу к лагранжевой координате. Ур. (3) в без- 
размерных переменных приведены в работах [4, 5]. 
В частности, весь процесс горения в модели опре- 
деляется следующими параметрами: 
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Здесь ноль вверху символов означает их стацио- 
нарные значения. Отрицательное / означает экзо- 
термическую реакцию на поверхности разложения 
топлива, положительное - эндотермическую. 


Результаты исследования и их анализ 


Как выяснилось в процессе решения ур. (2, 3), 
частота реализующихся в системе колебаний имеет 
порядок обратной величины т.н. аппаратурной 
константы X- 


_ К (»°) 2 
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представляющей собой отношение % = І ѵ /і„ где 
и^Ѵ с /(Я я Т р А с Р т ), І=к с /(и°) 2 - характерные времена 
истечения газа из сопла и тепловых процессов в 
твердой фазе топлива. Стационарная скорость го- 
рения м° определяется рабочим давлением р° в ка- 
мере,/- ІО 7 Па. 

Существует критическое значение %=%, ниже 
которого стационарный режим работы двигателя 
невозможен [1]. Как удалось установить, это кри- 
тическое значение / сильно зависит от параме- 
тров, определяющих механизм горения. Причем, 
уменьшение его значения означает появление все 
более высокочастотных колебаний с началом поте- 
ри устойчивости (рис. 1-4). Для всех рисунков 
ст=0,01; Іх=1,0; у=1,4. На рис. 1 /«2,4. Феномено- 
логические коэффициенты к, г рассчитывались по 
формулам, приведенным в [4, 5]. 

За границей устойчивости ^</«2,4 сначала на- 
ступает автоколебательное горение (рис. 1, а). По 
мере дальнейшего продвижения вглубь области не- 
устойчивости, т. е. с уменьшением %, колебания 
становятся более жесткими, с глубокими падения- 
ми скорости горения. В дальнейшем выхода на ав- 
токолебание не происходит: амплитуда скорости 
горения становится настолько большой, что реаги- 
рующая среда переходит к другому устойчивому со- 
стоянию - отсутствию горения. Проще говоря, на- 
ступает самопроизвольное погасание (рис. 1,6). 



Рис. 1. Потеря устойчивости, выход на автоколебательное 
горение (а, % =7,0 . ) и погасание (б, /=0,08): Ѳ 0 =О,3 ; 
І=-0,3 ; р=2, 0 ; г=0, 09; к=1, 04 

Как известно [6], горящее топливо как само- 
стоятельная система при неизменных внешних 
условиях имеет свой внутренний механизм, регу- 
лирующий его устойчивое состояние. Если этот 
механизм не сбалансирован, то наступает неустой- 
чивое горение, приводящее или к автоколебатель- 
ному горению, или к погасанию [4]. Удаленность от 
границы устойчивости (назовем ее С 0 -границей) 
может быть охарактеризована двумя параметрами к 
и г. Замечено, что величина / тем больше, чем 
дальше располагается состояние горения от 6і - гра- 
ницы. Например, стоит только немного отодви- 
нуться (по сравнению с рис. 1) от этой границы, 
наблюдается быстрое снижение / (рис. 2). Здесь 
уже /«0,63. 

Дальнейшее удаление от О/границы приводит 
к еще большему снижению /. Такое уменьшение 
сопровождается появлением новых качественных 
свойств: относительное небольшое изменение (в 
сторону уменьшения) % в неустойчивой области 
сразу же приводит к погасанию (рис. 3). Т. е. 
область существования по параметру % автоколеба- 
тельного режима горения сужается. Периодиче- 
ский режим горения сосредоточен в очень узком 
интервале изменения %, например, при ^=8- КТ 3 го- 
рение еще устойчиво. 

Варьирование параметром % показало возмож- 
ность реализации нерегулярного режима горения, 
который может быть охарактеризован как динамиче- 
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ский хаос (странный аттрактор). Действительно, по 
мере уменьшения аппаратурной константы растет 
амплитуда автоколебаний и все сильнее проявляется 
нелинейность системы. Если еще уменьшать %, то 
происходит удвоение периода колебаний. Над увели- 
чением периодов в дальнейшем проследить трудно: 
малое изменение % приводит к нерегулярному режи- 
му горения, подобно приведенному на рис. 4. Таким 
образом, переход к хаотическому режиму горения со- 
ответствует сценарию Помо-Манневила [7]. При 
больших значениях х хаотического колебания не 
удалось обнаружить. Хотя, конечно, это не говорит о 
его невозможности появления. Возможно, варьиро- 
вание параметрами произведено недостаточно. 

1,8 

1,6 

1,4 

1,2 

1,0 

0,8 

4 6 8 10 12 г 14 

1,6 
В 
1,4 

1,2 

1,0 

0,8 

0,6 

4 6 8 10 12 г 14 

Рис. 2. Потеря устойчивости и выход на автоколебательное 
горение: %=0,2; Ѳ 0 =О,3; І=~0,23; ц=2,2; т=0,13; 
к=1,17. Обозначены: Е,, - положение фронта пламени 
в газе ; Ѳ;, У 5 ~ температура и концентрация на по- 
верхности разложения 

Для надежной идентификации странного ат- 
трактора следовало бы произвести расчет энтропии 
Колмогорова-Синая или спектра реализующихся 
колебаний. Но это отдельная и довольно сложная 
задача. Но грубо наличие динамического хаоса 
можно определить визуально, представив колеба- 
тельный режим горения в фазовых переменных, 
например, в плоскости В - сІВ/Вт. Если зависи- 
мость В(т) периодическая функция, то на фазовой 
плоскости ей будет соответствовать замкнутая кри- 
вая. Если же периодичности нет, то на фазовой 
плоскости траектория будет «заметать» область [7]. 
Соответствующее построение приведено на рис. 5. 






0 2 4 6 8 10 


Рис. 4. Нерегулярные колебания скорости горения (а) и ос- 
новных параметров границ разделов фаз (б) при эн- 
дотермической реакции разложения топлива: 
Ѳ 0 =О,43 ; 1=0,23 ; р=3,0; к=0,70; т=0, 12; х=4 ■ Юг 3 
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Рис. 5. Фазовая кривая режима горения, приведенного на рис. 4 

Такому «визуальному» анализу аналогичен по 
своей сути подход, основанный на построении ото- 
бражения Пуанкаре: периодической кривой на 
плоскости Пуанкаре отвечает изолированный на- 
бор точек, а странному аттрактору - множество то- 
чек, которое может быть охарактеризовано своей 
размерностью, обычно дробным числом. 

Потеря устойчивости не всегда приводит к ав- 
токолебательному горению. Существуют такие па- 
раметры, например, Ѳ (] =0,6; /=-0,23; д=2,41; 
Д=0,06; /3=0,02; ^=3,8-10 8 ; а=0,01; Бе=1,0; у= 1,4; 
г= 0,14; к= 0,95, когда за границей устойчивости х 
возрастающие возмущения приводят к погасанию. 

Говорить о какой-либо конкретной частоте коле- 
баний на рис. 4, 5 не приходится. Но интервал (без- 
размерного) времени между двумя максимумами 
колебаний составляет около 0,07, т. е. сравним с без- 
размерным временем релаксации ст процессов в га- 
зовой фазе. Если этот промежуток времени принять 
за период колебаний, то при давлении около 10 7 Па 
и в размерных единицах их частота составляет при- 
мерно 10 3 Гц, т. е. относится к области высоких ча- 
стот. Для сравнения заметим, параметры на рис. 1, 2 
имеют частоту колебаний примерно 40 и 100 Гц. 

Заключение и основные выводы 

Моделирование нестационарного горения 
твердых ракетных топлив в камерах ракетных дви- 
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гателей показало возможность проведения описа- 
ния рассматриваемого явления на основе ур. (2, 3). 
Возникновение потери устойчивости при умень- 
шении аппаратурной константы, реализация вслед 
за этим автоколебательного и нерегулярного режи- 
ма горения или же погасания горения наблюдалось 
в многочисленных экспериментах, которые хоро- 
шо отражены в обзорной литературе [1, 8, 9 и др.]. 
Показанная в работе единая акустическая природа 
низкочастотной и высокочастотной неустойчиво- 
стей согласуется со сложившимися современными 
представлениями [1]. 

Возникновение акустической неустойчивости 
может приводить (если нет погасания) к росту или 
небольшому падению среднего давления в камере. 
В проведенных расчетах среднее давление всегда 
немного падает, несмотря на рост средней скорости 
горения. Это происходит из-за запаздывания давле- 
ния относительно скорости горения (рис. 1, а), т. к. 
изменение скорости горения зависит от давления 
опосредованно через изменение температуры. От- 
сюда следует вывод, что среднее давление в камере 
будет повышаться при снижении сдвига фаз между 
давлением и скоростью горения. Показатель ѵ 0 в 
скорости горения в приведенных выше расчетах по- 
лагался равным нулю. В противном случае принци- 
пиальных изменений в физической картине разви- 
тия неустойчивости не наблюдается, но ярче про- 
является нелинейный характер возникающих коле- 
баний. Это обычно приводит к нерегулярным изме- 
нениям во времени физико-химических параме- 
тров, характеризующих горение. 

Более подробное сравнение с имеющимися экс- 
периментальными данными не представляется воз- 
можным из-за отсутствия в литературе сведений по 
кинетике химических реакций, протекающих в ка- 
мерах ракетных двигателей. В настоящей работе 
при проведении расчетов ориентиром служили ко- 
эффициенты к а г, числовые значения которых ме- 
няются для большинства топлив в относительно не- 
больших интервалах [1, 6]: г=0,05...0,3; к= 0,8. ..1,5. 
Соответственно, параметры 0 О , /, д, Д Д, К () брались 
так, чтобы коэффициенты к и г попадали в указан- 
ные интервалы или же были близки к ним. 
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